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Abstract 

We perform a quenched lattice calculation of the first moment of twist-two generalized parton 
distribution functions of the proton, and assess the total quark (spin and orbital angular momen- 
tum) contribution to the spin of the proton. 
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Generalized parton distributions [l| (GPDs) provide a deeper understanding of the in- 
ternal structure of hadrons in terms of quarks and gluons. While ordinary parton dis- 
tributions measure the probability of finding a parton with fractional momentum 
X in the hadron, GPDs describe the coherence of two different hadron wave functions 
+ S,/2)ip{x — ^/2), one where the parton carries fractional momentum x + ^/2 and 
one where this fraction is x — ^/2, from which information about parton-parton correlation 
functions can be deduced. As a consequence, GPDs depend on the momentum transfer 
between the initial and final hadron, which provides further information on the transverse 
location of quarks and gluons j^. Spatial images of hadrons can thus be obtained, where the 
resolution is determined by the virtuality of the incoming photon. Last, but not least, 
GPDs allow us to isolate the contribution of the quark orbital angular momentum to the 
spin of hadrons. Lattice QCD is the only known method that is able to compute moments 
of GPDs from first principles. 

We will restrict ourselves to the GPDs Hq and Eq of the nucleon, where q = u,d, - ■ ■ 
denotes the flavor of the struck quark. We will not consider the gluon sector here. The 
lowest, zeroth moments of Hg and Eg are given by the Dirac and Pauli form factors: 



Both form factors have been computed on the lattice in a similar calculation 3] to the 
present one and found to be well described by a dipole ansatz 



for sufficiently small (and accessible) momenta, with dipole masses Mi^2 of the order of the 
p, u mass, when extrapolated to the physical pion mass. 

The first moments of Hg and Eg are of the form 




(1) 



(2) 



(3) 
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TABLE I: Parameters of the dipole fit. In the bottom row we give the parameters extrapolated to 
the physical pion mass. 





M [GeV] 


^"(0) 


Bm 


cm 


Ai{0) 


Bm 


cm 


0.1324 


1.69(05) 


0.419(07) 


0.344(028) 


-0.084(26) 


0.188(04) 


-0.281(20) 


-0.071(15) 


0.1333 


1.58(06) 


0.415(10) 


0.334(044) 


-0.101(35) 


0.176(05) 


-0.260(29) 


-0.073(19) 


0.1342 


1.41(10) 


0.404(19) 


0.357(117) 


-0.117(70) 


0.158(10) 


-0.265(80) 


-0.067(35) 




1.11(20) 


0.400(22) 


0.334(113) 


-0.134(81) 


0.147(11) 


-0.232(77) 


-0.071(42) 



where A\{A'^), i?|(A^) and C|(A^) are generahzed form factors, which are given by the 
nucleon matrix elements of the energy-momentum tensor (EMT) : 



{P'\0l^.^\p) = -{p'\q-i{,D^.}q\p) 



+ C|(A2 



2m 



-u{p')A''aa{^,Pu}u{p) 



(6) 



rriN 



N 



u{p')u{p) A{^A^} . 



Here rriN denotes the nucleon mass, p = 2 (p' + P); ^ = p' ~ and curly brackets refer to 
symmetrization of indices and subtraction of traces. The EMT has twist two and spin two. 
It is assumed to be renormalized at the scale fi, which makes A^IA"^), i?2(A^) and C|(A^) 
scale and scheme dependent. For the classification of states of definite J^'-' contributing to 
(jni) in the t-channel see The so-called skewedness parameter ^ is defined by ,^ = —n ■ A, 
where n is a light-like vector with n ■ p = 1, and bounded by |^| < 2^A^/{A'^-Am%). In 
the forward limit, A^ 0, we have 



^2(0) = {xq) = dxx (g|(x) + qi{x)) 



(7) 



where g-f(|)(x) are the usual quark distributions with spin parallel (antiparallel) to the spin 
of the nucleon. Furthermore, one derives 



l{Am+Bm) = jg, 



(8) 



where Jq is the angular momentum of the q quark, and J = Jq is the total angular 
momentum of the nucleon carried by the quarks. The angular momentum decomposes, in a 
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gauge invariant way, into two pieces: 

Jg=Lg + Sq, (9) 

where Lg is the orbital angular momentum and 

Sg = ^Aq^^l^ dx {q^{x) - q^ix)) (10) 

is the spin of the quark. We know Aq from separate calculations 0, Ql , so that Lq can be 
computed from (jS)). 

In this Letter we perform a quenched lattice calculation of the generalized form factors 
AKA"^), i?|(A^) and C|(A^). The quenched approximation neglects fluctuations of vir- 
tual quark-antiquark pairs from the Dirac sea. The non-forward matrix elements (EI) are 
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computed from ratios of three- and two-point functions following p|. Further details are 
given in |8|. To keep cut-off effects small, we use non-perturbatively 0{a) improved Wilson 
fermions. The calculation is done on 16^ 32 lattices at f3 = 6.0 and for three different hop- 
ping parameters, n = 0.1324, 0.1333 and 0.1342, which allows us to extrapolate our results 
to the chiral limit. Using ro = 0.5 fm to set the scale, which results in the inverse lattice 
spacing 1/a = 2.12 GeV, the corresponding pion masses are 1070, 870 and 640 MeV. If we 
use the nucleon mass extrapolated to the chiral limit to set the scale, the pion masses are 
930, 760 and 550 MeV, and 1/a = 1.84 GeV. The corresponding nucleon masses and the 

n 

choice of nucleon momenta p,p' can be inferred from [3]. For the EMT we consider two sets 
of (euclidean) operators: 

-l={0,,, + 0,,,), l<fl<U<4: (11) 

and 

^(011 + 022-033 -(^44), 

(12) 

-^(033-044), ^(011-022). 

Each set transforms irreducibly under the hypercubic group. The operators PHI and (fT^ are 
renormalized multiplicatively, 0{fi) = Z{afi) 0{a), with renormalization constants |^ 
and Z„26! respectively. The renormalization constants are computed non-perturbatively pf 
following [lo|. We obtain ^^^(2 GeV) = 1.10 and Z^(2 GeV) = 1.09. The following 



results refer to the MS scheme at the renormalization scale /x = 2 GeV. 
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FIG. 1: The generalized form factors A2, B2 and at k = 0.1333, together with the dipole fit 

and the extrapolated values at = (□). 
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FIG. 2: The generalized form factors Ag, i?2 and at k = 0.1333, together with the dipole 

and the extrapolated values at = (□). 



6 



In Figs. 1 and 2 we show the generahzed form factors Al^iA.'^), i?2(A^), C2(A^) and 
y42(A^), i?2(A^), C2(A^) of the proton for k = 0.1333. Data points with larger errors are 
not shown here but are included in the fit. The corresponding form factors of the neutron 
are obtained by interchanging u and d. Similarly good results are found for n = 0.1324 and 
0.1342. The generalized form factors can be well described by the dipole ansatz 

Al{A') = Al{0)/{l-AyMy, (13) 

and similarly for B2 and C|. Fits of y42(A^) and A^IA"^) give the same dipole mass M 
within errors. The dipole masses obtained from separate fits of i?2(A^), i?2(A^), C2(A^) 
and 6*2 (A^) are found to be consistent with that value. We therefore have decided to fit our 
data by a common dipole mass M. Our data do not favor a monopole behavior. The results 
of the fits are shown in Table I. For a reliable extrapolation to A^ = we find it important 
to cover a wide enough range of A^ values. This may be the reason why our dipole masses 
turn out to be systematically larger than those found in a previous calculation jll| . 

In Fig. 3 we show the dipole mass M as a function of the pion mass. The mass values 



appear 
factors 



;o lie on a straight line, as was observed already in the case of the nucleon form 



a 



A linear extrapolation in to the physical pion mass gives M = 1.1(2) GeV. 
This value is close to the physical masses of the /2, 02 mesons, which supports the hypothesis 
of tensor meson dominance. A quadratic extrapolation in m,r leads to M = 1.3(1) GeV. 
The form factor data Al{0), -62(0) and C|(0) show little variation with the quark mass 
and are extrapolated quadratically in m^r to the physical pion mass. The results are shown 
in the bottom row of Table I. It should be stressed that all quantities refer (at best) to 
valence quark distributions, because sea quark effects have been neglected. In unquenched 
simulations there are also quark-line disconnected contributions. For an estimate see jll| . 

If the dipole behavior Q, ()13p continues to hold for the higher moments as well, and if 
we assume that the dipole masses continue to grow in a Regge-like fashion, we would obtain 

l\xx-H,{x,0,A') = (x^)/(l- AVM„V)', (14) 

with = const. + l/a', a' being the slope of the Regge trajectory. This would mean that 
with increasing momentum transfer |A^| the lower moments of iJg(x,0, A^) are suppressed 
more than the higher ones, so that the observed peak in Hg{x, 0, 0) = qi^{x) + qi{x) around 
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FIG. 3: The dipole mass M as a function of ttt-tt, together with a linear extrapolation to the physical 
pion mass (□). 

X ~ 0.2 is shifted towards the higher values of x. As a result, the dependence cannot 
be factorized in a simple way, as is sometimes assumed. Knowing (xp, we can reconstruct 
Hq{x,0, A"^) from (fT^ by inverse Mellin transform. The ^ dependence of both Hg and 
Eg appears to be rather weak, based on our knowledge of the first two moments, and in 
the isovector channel (corresponding to proton-neutron or u-d matrix elements) it largely 
cancels out. 



In Fig. 4 we show the total angular momentum J = 3^ + Jd- The dependence on the pion 
mass is rather flat, as expected The errors are due to the relatively large statistical 

errors of B2 and B2 and the fact that B2 and B2 cancel to a large extent. In Table II we give 
our results for J, and separately for Jg and Sq, extrapolated quadratically (linearly in m^) 
to the physical pion mass. The numbers for Sq refer to our latest results [9] , computed from 
the non-perturbatively improved axial vector current with non-perturbative renormalization 
factors. It turns out that the total angular momentum J carried by the quarks amounts to 
~ 70% of the spin of the (quenched) proton, leaving a contribution of ~ 30% for the gluons. 
The major contribution is given by the u quark, while the contribution of the d quark is 
found to be negligible, which hints at strong pairing effects. Our result for J is somewhat 



smaller than that of 



. We are able to compute L„ now. The total orbital angular 
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FIG. 4: The total angular momentum J, together with a quadratic extrapolation to the physical 
pion mass (□). 



TABLE II: The total angular momentum and its individual contributions, extrapolated to the 
physical pion mass. 



J 




Jd 




Sd 


0.33(7) 


0.37(6) 


-0.04(4) 


0.42(1) 


-0.12(1) 



0.4 



momentum of the quarks turns out to be consistent with zero: 

L = L„ + Lrf = 0.03(7). (15) 

This indicates that (at virtuality fi = 2 GeV) the parton's transverse momentum in the 
(quenched) proton is small. A similar conclusion can be drawn from our earlier finding 
of a small twist-three contribution ^2 to the second moment of the polarized structure 
function g2. 

The generalized form factors C|(A^) contribute to the beam charge asymmetry of deeply 
virtual Compton scattering. We obtain a rather small value: 6*2(0) +6*2(0) = —0.2(1). This 
result is to be compared with the value —0.8 obtained in the chiral quark-soliton model at 
~ 0.6 GeV |l^. For a discussion see also 

As far as one can compare, quenched and unquenched results agree surprisingly well, and 
we do not expect to find significant differences here either. For a recent study of quenching 



artifacts, as well as cut-off effects, see jl7[ . 
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